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ABSTRACT 

Both astronomical observations of the interaction of Type II snpernova rem¬ 
nants (SNR) with dense interstellar clonds as well as cosmochemical stndies of the 
abnndances of danghter prodncts of short-lived radioisotopes (SLRIs) formed by 
snpernova nncleosynthesis snpport the hypothesis that the Solar Systems SLRIs 
may have been derived from a snpernova. This paper continnes a series devoted 
to examining whether snch a shock wave conld have triggered the dynamical 
collapse of a dense, presolar clond core and simnltaneonsly injected snfficient 
abnndances of SLRIs to explain the cosmochemical evidence. Here we examine 
the effects of shock waves striking clonds whose spin axes are oriented perpen¬ 
dicular, rather than parallel, to the direction of propagation of the shock front. 

The models start with 2.2 Mq cloud cores and shock speeds of 20 or 40 km s“^. 
Central protostars and protoplanetary disks form in all models, though with disk 
spin axes aligned somewhat randomly. The disks derive most of their angular 
momentum not from the initial cloud rotation, but from the Rayleigh-Taylor 
hngers that also inject shock wave SLRIs. Injection efficiencies, fi, the fraction 
of the incident shock wave material injected into the collapsing cloud core, are 
0.04 - 0.1 in these models, similar to when the rotation axis is parallel to the 
shock propagation direction. Evidently altering the rotation axis orientation has 
only a minor effect on the outcome, strengthening the case for this scenario as 
an explanation for the Solar Systems SLRIs. 

Subject headings: hydrodynamics — instabilities — ISM: clouds — ISM: super¬ 
nova remnants — planets and satellites: formation — protoplanetary disks — 
stars: formation 
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1. Introduction 

Radio and sub-millimeter observations of the Type II supernova (SNe) remnant (SNR) 
W44 and its interaction with the W44 giant molecular cloud complex imply the presence 
of dense clumps of molecular gas that have been struck and compressed to sizes much less 
than 0.3 pc by the SNR shock front, with shock speeds of 20 to 30 km s“^ (Sashida et 
ah 2013; Reach et ah 2005). HST images of the Cygnus Loop Type II SNR reveal shock 
front thicknesses of order 2 x 10“"^ pc (Blair et ah 1999). Detailed two dimensional (2D) 
axisymmetric (Boss & Reiser 2013) and fully three dimensional (3D) (Boss & Reiser 2012, 
2014; Li et ah 2014) hydrodynamical models of shock waves with these speeds and thicknesses 
have demonstrated the viability of triggered collapse of dense cloud cores with simultaneous 
injection of shock front material into the resulting collapsing protostar and protoplanetary 
disk. Given the problems with shocks waves associated with AGB and WR star winds (Boss 
& Reiser 2013), a SNR shock appears to be the leading contender for achieving simultaneous 
triggering and injection of core-collapse supernova-derived short-lived radioisotopes (SLRIs, 
e.g., ®°Fe and ^®A1) into the presolar cloud and the resulting solar nebula (Gameron & Truran 
1977; Boss 1995). 

Gosmochemical evidence for live ®'^Fe during the formation of refractory inclusions and 
chondrites has waxed and waned in recent years. Data on ferromagnesian chondrules from 
two ordinary chondrites (OG) implied initial ®°Fe/®®Fe ratios of ~ 5 — 10 x 10“^ (Tachibana 
et al. 2006), whereas bulk sample data from a wide range of meteorites suggested an initial 
ratio of only ~ 1.15 x 10“® (Tang & Dauphas 2012). A combined study of ®°Fe and ^®A1 in 
chondrules from unequilibrated OG (UOG) implied an initial ®°Fe/®®Fe ratio of ~ 7 x 10“^ 
(Mishra & Goswami 2014) and supported a SNe as the source of the SLRIs. Analyses of other 
chondrules from UOG yielded ratios in the range of ~ 2 — 8 x 10“^ (Mishra & Ghaussidon 
2014). While there is as yet no explanation for the discrepancy between bulk samples and 
chondrule fragments (Mishra & Ghaussidon 2014), a SNe remains as a plausible source for the 
SLRIs. Indeed, evidence for possible live ^^^Gs in GAIs suggests its origin in a Type II SNe (or 
its WR progenitor) close to the presolar cloud (Bermingham et al. 2014). Furthermore, the 
evidence for live ^°Be in FUN-GAIs appears consistent with ^'^Be formation in the presolar 
molecular cloud by galactic cosmic rays (GGRs) emitted by a SNR that thereafter triggered 
the collapse of the presolar cloud core (Tatischeff et al. 2014). Statistical analysis of the SLRI 
enrichments expected for disks and cloud cores in clusters containing massive stars showed 
that cloud cores receive a larger dose on average than disks (Adams et al. 2014), favoring the 
presolar cloud scenario (Gameron & Truran 1977; Boss 1995) over the alternative mechanism 
of injection into a previously formed solar nebula (Ouellette et al. 2007, 2010). 

Given this strong support for the triggering and injection hypothesis for the origin of 
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the SLRIs, we continue to calculate increasingly detailed models of the hydrodynamics of 
the shock-cloud interactions, in order to further understand the physics of the process. Our 
previous work in both 2D (Boss & Keiser 2013) and 3D (Boss & Keiser 2014) considered 
rotating presolar clouds, capable of collapsing to form protostellar disks. Because of sym¬ 
metry constraints in 2D, these models assumed that the cloud rotation axis was parallel to 
the direction of shock wave propagation. For ease of comparison, our first 3D rotating cloud 
models (Boss & Keiser 2014) made this same assumption. Here we present the results of a 
set of 3D models identical to previous rotating cloud models, but with the rotation axis per¬ 
pendicular to the shock front direction. We seek to discover what effect, if any, the rotational 
axis orientation might have on the SNe trigger and injection hypothesis. 


2. Numerical Hydrodynamics Code 

The numerical models were calculated with the FLASH2.5 adaptive mesh refinement 
(AMR) hydrodynamics code (Fryxell et al. 2000) in exactly the same manner as in our 
previous 3D models (Boss & Keiser 2012, 2014), except for changing the initial direction of 
the target cloud rotation axis to be perpendicular to the shock propagation direction. We 
give here a brief summary of the basic details of our implementation of the FLASH2.5 AMR 
code; further details about our implementation may be found in Boss et al. (2010). 

We employed FLASH2.5 in Cartesian coordinates (x, y, z). FLASH uses a block-structured 
adaptive grid approach, using the PARAMESH package, to enhance the grid resolution in 
regions with strong gradients of chosen variables. Initially the numerical grid has six grid 
blocks along the x and i axes and nine along the y axis. The —y axis is the direction of 
propagation of the shock wave, while the rotation axis of the target cloud is the x axis. The 
grid blocks each consist of 8^ grid points with four levels of grid rehnement initially. Rehned 
grids have one half the grid spacings of their parent grids. We use both the density and 
the color held (dehned below) to determine when grid rehnement is necessary to follow the 
shocked regions. While Boss et al. (2010) found that adequate resolution was obtained with 
a maximum of hve or six levels of grid rehnement, here the maximum number of levels of grid 
rehnement permitted was increased in some models to seven levels during the evolutions, 
when the memory allocation of the hash cluster nodes used for the calculations permitted 
the use of yet another level of grid rehnement. Additional levels of grid rehnement (beyond 
four) are necessary in order to attain the highest possible resolution of the dynamically col¬ 
lapsing regions that occur late in the evolutions. With seven levels, the smallest grid spacing 
reached is then ~ 4 x 10“^ pc ~ 9 AU. 

As in our previous models (Boss et al. 2010; Boss & Keiser 2012, 2013, 2014), we used a 



4 


color field to follow the evolution of shock front gas and dust. This allows us to quantify the 
efficiency with which the shock front matter is injected into the target molecular cloud core. 
The target molecular cloud core and the surrounding, ambient low density gas are initially 
isothermal at 10 K. The shock front and post-shock gas are initially isothermal at 1000 
K. This initial shock temperature follows from the C-shock models of Kaufman & Neufeld 
(1996), who found peak shock temperatures of ~ lO^K for shock speeds in the range of 5-45 
km s“^ (see Boss et al. 2010). We then follow the compressional heating as well as cooling 
of the interacting shock-cloud system by molecular species. In particular, we assume cooling 
at a rate consistent with cooling by optically thin H 2 O, CO, and H 2 molecules in order to 
follow the thermodynamics of the shock-cloud interaction. As before, we used Neufeld & 
Kaufman’s (1993) radiative cooling rate of A ^ 9 x 10^^(T/100)p^ erg cm“^ s“^, where T is 
the gas-dust temperature in K and p is the gas density in g cm“^. This rate is appropriate 
for cooling caused by rotational and vibrational transitions of optically thin, warm molecular 
gas (Boss et al. 2010). Note that we do not try to follow the collapsing protostellar objects 
and disks deep into the nonisothermal collapse regime, after the collapsing central regions 
become optically thick at densities above ~ 10“^^ g cm“^; we leave this additional physical 
complication for onr futnre models. 


3. Initial Conditions 

Table 1 lists the variations in the initial conditions for the six new 3D models. These 
models have all been previously calculated with the rotational axis aligned with the shock 
propagation direction, allowing a direct comparison of the resnlts in each case (Boss & 
Keiser 2014). All models start with 2.2 Mq molecular cloud cores with radii of 0.053 pc and 
Bonnor-Ebert radial density prohles. The clouds have solid body rotation abont the x axis 
with an angnlar frequency of either Dc = 10“^^ or 10“^^ rad s“^. These angular frequencies 
are typical of observed dense clond cores (e.g., Harsono et al. 2014). 

The shock waves propagate toward the —y direction with shock speeds of Vg = 20 or 40 
km s“^. The shock widths are 3 x 10“^ pc for all the models, i.e., the shock widths Wg = 0.1 
in units of the standard shock width, 0.0030 pc (Boss et al. 2010). The initial shock densities 
{ps) are 7.2 x 10“^® g cm“^ or 1.44 x 10“^^ g cm“^, i.e., 200 or 400 times the reference shock 
density of 3.6 x 10“^° g cm“^. Details abont why these particular shock wave parameters 
are appropriate for a supernova shock that has travelled several pc through the interstellar 
medinm prior to smacking the target clonds may be fonnd in Boss & Keiser (2012, 2013, 
2014). In particular. Boss et al. (2010) studied the effects of shock speeds ranging from 1 
km s“^ to 100 km s“^, and fonnd that only shock speeds in the range of 5 km s“^ to 70 km 
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s“^ were able to trigger collapse and injection into the target cloud assumed here. Boss & 
Keiser (2013) studied the effects of different shock widths {wg = 0.1, 1.0, and 10.0). Our use 
of tCs = 0.1 is consistent with observations of the Cygnus Loop (Blair et al. 1999) and other 
SNRs (see Boss & Keiser 2014). 


4. Results 

As before (Boss & Keiser 2012, 2013, 2014), the shocks strike the top edge of the 
target clouds and compress the central regions into sustained self-gravitational collapse, 
while accelerating the clouds downstream in the shock direction. The molecular line cooling 
keeps the temperatures of the shocked region below 1000 K, while the un-shocked regions 
remain isothermal at 10 K. Rayleigh-Taylor (R-T) hngers form at the shock-cloud interface, 
and these R-T hngers are directly responsible for injecting supernova-derived material into 
the collapsing cloud core. The calculations are terminated once the maximum densities 
exceed ~ 10“^^ g cm“^, as at those densities the collapsing regions have become optically 
thick and will no longer be able to cool at the rate appropriate for optically thin regions. 

Figures 1 through 6 show the basic results for the six models dehned in Table 1. It is 
evident that central protostars and edge-on disks form in all six models, though with overall 
disk structures considerably different than in models where the shock direction is parallel 
to the cloud rotation axes (Boss & Keiser 2014). In the latter case, the edge-on disks are 
always aligned with their innermost midplanes perpendicular to the shock direction, i.e., 
with their rotational axes aligned with the shock direction (cf.. Figure 4 of Boss & Keiser 
2014). This is a direct consequence of their initial rotational axes pointing downward, in the 
sense of Figures 1 through 6. The outermost disks, however, in the aligned axes models do 
warp upwards or downwards compared to the innermost disk midplane, beyond distances 
of ~ 150 AU for ric = 10“^"^ rad s“^ and beyond ~ 500 AU for hlc = 10“^^ rad s“^. In 
comparison, the edge-on disks seen in Figures 1 through 6 in some cases have innermost disk 
midplanes tilted compared to the perpendicular to the shock direction (e.g.. Figures 3, 4, 5), 
while the remainder are closer to being perpendicular (Figure 1, 2, and 6). In all six models, 
however, the outermost disks warp away from the midplane orientations. 

What is most remarkable about the disk orientations in Figures 1 through 6 is the fact 
that in no case are the disk midplanes seen to be aligned with the direction of the shock 
propagation, i.e., the —y axis, as might be expected considering that the initial sense of 
rotation of the target clouds is around an axis aligned with the +x axis. This is quite 
different from the results for the Boss & Keiser (2014) shock-aligned axes models, where the 
disk midplanes are perpendicular to the initial cloud spin axes, implying that the disks have 
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derived their angular momentum from the initial target clouds. Given the fact that the new 
disks are warped and otherwise highly distorted (see especially Figure 5), it is unclear how 
to assign a single, specihc direction of an overall disk rotational axis, as there evidently is no 
single rotational axis for these models. Considering the situation evident in Figures 
1 through 6, Table 1 lists the variety of spin directions ascertained for different 
portions of each disk. 

The reason for this puzzling difference in the results can be discovered hrst by examining 
the spin directions of the material in the disks. Figure 7 shows the component of the velocity 
held for model 40-200-prl3, plotted in the same sense as in Figure 1 for this model. Figure 
7 shows that the axis of rotation of the innermost disk is along the +y axis, in the opposite 
direction of the shock propagation, and perpendicular to the initial target cloud rotation axis. 
Figure 8 shows the same type of plot for model 20-400-prl4, where the axis of rotation is now 
aligned with the —y axis, the opposite direction of model 40-200-prl3, again perpendicular 
to the initial target cloud rotation axis. Table 1 summarizes the approximate disk spin 
axis directions for all six models, showing no obvious pattern for the overall orientations. 
Evidently the disks seen in Figures 1 through 6 did not derive the majority of their angular 
momentum from the assumed initial cloud angular momentum. 

Further investigation of Figures 7 and 8 reveals the source of the bulk of the angular 
momentum contained in the disks. Figure 7 shows that the largest variations in the velocity 
held are found in the R-T hngers in the shock front. In fact, in Figure 7, a clear stream of 
positive Vz gas can be traced from the R-T hngers down to the innermost disk midplane. 
Similarly, in Figure 8 a large held of negative Vz gas connects the disk midplane with the 
R-T hngers. This explanation is made even more robust in Figures 9 and 10, where the disk 
formed in model 40-200-prl3 is seen face-on, with a strong trailing one-arm spiral indicative 
of the disk spinning with its rotational axis aligned with the +y direction. Figure 10 plots 
the Vx velocity component for this model’s midplane, showing the same sense of rotation as 
previously inferred. A plot of Vz shows the same agreement. Again, Figure 10 shows that the 
greatest variations in Vx originate in the R-T hngers, and in fact streams of similar velocity 
gas in the disk can be traced back to R-T hngers with the same velocity sign. Evidently 
the R-T hngers provide the bulk of the angular momentum for disks formed when the initial 
rotational axes are not aligned with the shock propagation direction. In general, one would 
not expect these two to be aligned, so it is important to understand how this could happen. 

To further investigate this odd result, it is important to note that when the initial 3D 
cloud is not rotating at all, i.e., when Dc = 0, no obvious disks form at all, in spite of 
R-T hngers being present (Boss & Reiser 2012). Instead, the collapsing protostar is roughly 
spherical (cf.. Figure 3 in Boss & Reiser 2012). This implies that it is the interaction of the 
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initial target cloud spin with the R-T hngers formed at the shock-cloud interface that leads 
to the R-T hngers being sheared and dehected away from motions that would otherwise be 
more purely directed downward, along the —y axis, he., the Vy momentum of the shock 
gets converted by the sheared R-T hngers into a complex pattern of and momentum, 
sufficient to provide the angular momentum of the disks seen in Figures 1 through 6. In 
support of this explanation, note that the Keplerian velocity of gas orbiting at 100 AU 
around a protostar with a mass of 0.5 Mq is ~ 2 km . Figure 7 shows that ~ 2 km is 
indeed a typical speed for the velocity component of the gas 100 AU out in the rotating 
edge-on disk in model 40-200-prl3, consistent with a disk in Keplerian rotation. Similarly, 
Figure 10 illustrates exactly the same situation for the x component of the gas velocity at 
100 AU from the central protostar in the face-on disk seen in Figure 9, for this same model. 
These speeds are the same as can be seen in the R-T hngers evident in Figures 7 and 10, from 
which the disk gas originated, and cement the case for the angular momentum sufficient to 
form rotating, Keplerian protostellar disks being derived from the R-T hngers, rather than 
from the initial cloud rotation. 

In this same context, it is important to note that the shocks start with considerably 
higher speeds along y (20 to 40 km s“^) compared to the initial target cloud rotational 
speeds (no more than ~ 0.1 km s“^ for hlc = 10“^^ rad s“^, and no more than ~ 0.01 km s“^ 
for hlc = 10“^^ rad s“^), so it is easy to see how even a small dehection of the shock front 
momentum carried by the R-T hngers can dominate the outcome for the disks that form. 
This phenomenon also explains the origin of the outer disk warps seen in the models with 
aligned rotational axes and shock directions (Boss & Keiser 2012): the addition of material 
to the disk with velocities inhuenced by the R-T hngers. 

This explanation appears to be consistent with the quite diherent result of Li et al. 
(2014), who showed one model (R2) where an initial cloud rotating similar to the present 
models produced a disk with a spin axis in the expected direction (cf.. Figure 2(d) in Li et al. 
2014). While there are many diherences between the present models and those of Li et al. 
(2014), e.g., their use of sink cells and our use of detailed thermodynamics, the fact that their 
R2 model started with a shock speed of 3 km s“^ implies that the shock-derived momentum 
was unable to overcome the initial cloud rotational motions. In addition, R-T hngers do not 
appear in their calculations, apparently because of insufficient spatial resolution, in spite of 
their using an AMR code, ASTROBEAR 2.0 (Li et al. 2014). However, Li et al. (2014) 
used only three levels of grid rehnement, leading to a smallest grid size of about 23 AU, 
whereas the present models with seven levels of rehnement had a smallest grid size of about 
9 AU, which could be an important factor. Combined with the other diherences between 
these two sets of models, the fact that Li et al. (2014) did not observe R-T hngers in any of 
their models ensures that, unlike the present models, where R-T hngers are essential to the 



injection process, R-T fingers could not dominate the angular momentum of the disk that 
formed in the Li et ah (2014) model R2. 

These calculations imply that besides injecting the SLRIs found in primitive meteorites, 
R-T fingers may have also delivered the angular momentum that allowed a protostellar disk 
to form around the nascent protosun, a disk that would evolve to become the protoplanetary 
disk that formed our Solar System. The fact that the disk spin axis orientations listed in 
Table 1 are somewhat random implies that protoplanetary disks formed by this mechanism 
would be expected to have spin axes with roughly random orientations, orientations that are 
determined by the detailed development of the R-T hngers that delivered the disks’ angular 
momenta. 

We now turn to a consideration of the injection efficiencies for these models. As before 
(Boss & Reiser 2012, 2013, 2014), the injection efficiency /* is defined as the fraction of the 
incident shock wave material that is injected into the collapsing cloud core. Table 1 lists 
the calculated values of fi, along with the mass (Me) of the dynamically collapsing regions 
(i.e., regions with p > 10“^® g cm“^) and the mass of the shock wave-derived material in this 
same region {Mg). The calculated values of fi compare favorably with those for otherwise 
identical models with aligned rotation (cf.. Table 1 in Boss & Reiser 2014). E.g., model 
40-200-prl3 has fi = 0.037, whereas the corresponding aligned rotation model from Boss & 
Reiser (2014) had fi = 0.032. This same minor increase in /* is seen for the other models 
in Table 1 with corresponding aligned axes models as well. Clearly having non-aligned axes 
does not make SLRI injection any harder. The increased values of fi can be attributed to 
the larger role of the R-T fingers in the non-aligned models in providing the bulk of the 
disk angular momentum, along with the SLRIs. Figures 11 and 12 show the color fields 
for models 40-200-prl3 and 20-400-prl4, making it clear that by this advanced phase of the 
evolutions, the R-T fingers have successfully injected shock wave material throughout the 
collapsing regions. 

While the new models exhibit the robustness of our previous estimates of fi, it must 
also be noted that the SNe shock waves must sweep up considerable amounts of intervening 
interstellar gas and dust in order to slow down to the shock speeds considered here and 
previously found to be the most favorable for injection (e.g.. Boss et ah 2010; Boss & Reiser 
2013). This implies a large amount of dilution of the pristine SLRIs formed by the SNe. 
As in Boss & Reiser (2012), one can dehne a factor to be the ratio of the shock mass 
originating in the SNe to the mass swept up in the intervening gas and dust. For the present 
shock speeds, a value of /3 ^ 0.01 is appropriate (Boss & Reiser 2012). The dilution factor 
D is the ratio of the amount of mass derived from the supernova to the amount of mass 
derived from the target cloud. For models with ps = 200 (Table 1), the amount of mass in 
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the incident shock wave is 0.30 Mq. Assuming a hnal mass for the protostar and disk of 
1.0 Mq, this results in H ~ 3 x 10“"^. This dilution factor is in reasonable agreement with 
estimates for SNe-derived SLRIs of H ~ 10“"^ to H ~ 10“^ (Takigawa et ah 2008). In spite 
of the ongoing uncertainty over the correct initial abundances of SLRIs in the solar nebula, 
as discussed in the Introduction, the dilution factor estimates of Takigawa et al. (2008) 
appear to be the best current values, as argued in Boss & Reiser (2014). Li et al. (2014) 
found insufficient SLRI injection in their models to explain the cosmochemical 
evidence. However, as noted above, Li et al. (2014) were not able to resolve 
the R-T fingers necessary for significant injection, and furthermore studied only 
shock speeds of 3 or 6 km s“^, at or below the lower end of the range (5 to 70 
km s“^) found necessary for simultaneous triggering and injection by Boss et al. 
( 2010 ). 

Finally, we note from Table 1 that the injection efficiencies fi can be seen to increase as 
the amount of mass Me in the collapsing region with p > 10“^® g cm“^ increases. Evidently 
some of the estimates of fi in Table 1 may be lower bounds, in that the values of fi are likely 
to increase when the calculation is taken farther in time and more mass joins the collapsing 
regions. In practice, the calculations are halted when the collapsing region becomes too dense 
for the isothermal approximation to be valid (i.e., p > 10“^^ g cm“^), which is also about 
when the collapsing regions reach the bottom edge of the computational domain. Both of 
these limitations are being addressed in a new set of models currently underway. 


5. Conclusions 

These models have demonstrated that the degree of alignment between a rotating molec¬ 
ular cloud core and the direction of propagation of a shock wave that strikes the cloud core, 
triggering its self-gravitational collapse, does not have a major effect on the efficiency of 
injection of shock-derived SLRIs into the collapsing cloud core. In fact, the injection effi¬ 
ciencies in the non-aligned models (~ 0.04 to 0.1) are somewhat higher than in the aligned 
models with identical shock wave parameters (cL, Boss & Reiser 2012, 2013, 2014). Hence 
from the point of view of the numerical calculations of the mechanism, the supernovae shock 
wave triggering and injection hypothesis remains as a plausible explanation (Cameron & 
Truran 1977; Boss 1995). 

Remarkably, these new models have also introduced a new feature of the shock wave 
triggering and injection mechanism; the R-T fingers responsible for SLRI injection can 
concomitantly result in the injection of enough momentum to largely determine the direction 
of the resulting disks’ spin axis orientations. In such cases, the R-T fingers may have been 
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responsible not only for the acquisition of the SLRIs inferred to have been present in the 
most primitive meteorites, but also for the very fact that a rotating protostellar disk was 
formed, a disk that eventually led to the formation of our planetary system. 

Our future 3D FLASH models will investigate the effects of the loss of molecular line 
cooling once the clouds become optically thick at densities well above ~ 10“^^ g cm“^. In 
this regime, the collapsing regions begin to heat above 10 K, but continue their collapse 
toward the formation of the first protostellar core, at pmax ~ 10“^° g cm“^ (e-g-. Boss & 
Yorke 1995). 


We thank the referee for improvements to the manuscript, and Michael Acierno for help 
with the flash cluster at DTM, where the calculations were performed. This research was 
supported in part by NASA Origins of Solar Systems grant NNX09AF62G. The software 
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Astrophysical Thermonuclear Flashes at the University of Chicago. 
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Table 1. Initial conditions and basic results for the models, with varied shock speeds 
in units of km s“^), shock gas densities {ps, in units of the standard shock density, 

3.6 X 10“^° g cm“^), target cloud perpendicular rotation rates (ffc, in rad s“^), mass in 
collapsing (p > 10“^® g cm“^) region {Me, in in Mq), shock wave masses in collapsing 
regions {Mg, in Mq), injection efficiencies (/*), and approximate spin directions of resulting 

disks. 


Model 

Vs 

Ps 



Me 

Ms 


spin direction 

40-200-prl3 

40 

200 

10“ 

-13 

0.45 

0.011 

0.037 

+y 

40-200-prl4 

40 

200 

10- 

-14 

0.49 

0.012 

0.040 

-y 

20-200-prl3 

20 

200 

10- 

-13 

0.95 

0.028 

0.093 

-y to +x to +z 

20-200-prl4 

20 

200 

10- 

-14 

1.04 

0.031 
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Fig. 1.— Model 40-200-prl3 log density cross-section in the z = —1.5 x 10^^ cm plane, where 
the maximnm density occnrs, showing an edge-on disk. Region shown is 1333 AU across at 
0.085 Myr. Shock propagates from top to bottom in each model. Direction of initial target 
clond spin axis is to the right. 
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Fig. 2.— Model 40-200-prl4 log density, plotted as in Fig. 1, bnt in the z = —0.5 x 10^® cm 
plane, at 0.079 Myr. 
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Fig. 3.— Model 20-200-prl3 log density, plotted as in Fig. 1, bnt in the z = —1.8 x 10^® cm 
plane, at 0.15 Myr. 
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Fig. 4.— Model 20-200-prl4 log density, plotted as in Fig. 1, bnt in the z = 0 plane, at 0.13 
Myr. 
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Fig. 5.— Model 20-400-prl3 log density, plotted as in Fig. 1, bnt in the z = —1.2 x 10^® cm 
plane, at 0.093 Myr. 
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Fig. 6.— Model 20-400-prl4 log density, plotted as in Fig. 1, bnt in the z = 0 plane, at 
0.090 Myr. 
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Fig. 7.— Model 40-200-prl3 2 ; velocity, plotted as in Fig. 1 at 0.085 Myr. 
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Fig. 8.— Model 20-400-prl4 2 ; velocity, plotted as in Fig. 6, at 0.090 Myr. 
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Fig. 9.— Model 40-200-prl3 log density in the disk midplane {y = 2.24 x 10^® cm) at 0.085 
Myr. 
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Fig. 10.— Model 40-200-prl3 x velocity in the disk midplane {y = 2.24 x 10^® cm) at 0.085 
Myr. 
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Fig. 11.— Model 40-200-prl3 color density, plotted as in Fig. 1 at 0.085 Myr. 
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Fig. 12.— Model 20-400-prl4 color density, plotted as in Fig. 6 at 0.090 Myr. 















